Tunneling theory of two interacting atoms in a trap 
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A theory for the tunnehng of one atom out of a trap containing two interacting cold atoms is 
developed. The quasiparticle wave function, dressed by the interaction with the companion atom 
in the trap, replaces the non-interacting orbital at resonance in the tunneling matrix element. 
The computed decay time for two ^Li atoms agrees with recent experimental results [G. Ziirn, F. 
Serwane, T. Lompe, A. N. Wenz, M. G. Ries, J. E. Bohn, and S. Jochim, arXi v:1111.2727| , unveiling 
the 'fermionization' of the wave function and a novel two-body effect. 

PACS numbers: 67.85.Lm, 03.65.Xp, 73.23.Hk, 03.65.Ge 




FIG. 1. Gonfinement potential V{x) vs x according to 
Ref.|6l The harmonic oscillator frequency ujq = 0.744 KHz-27r 
is obtained from the spacing between the two bound states of 
V{x) computed through the WKB approximation (solid thin 
lines); £ is the resonating tunneling energy. The length unit 
is feo = {h/muJoY^'^ . The circles with arrows schematize the 
atoms in the W\) configuration. 



Recently few cold atoms have been confined in tight 
single optical traps with control of their number and 
quantum states [P, '2^. This capability is exciting as it 
bridges different fields such as mesoscopic and nuclear 
physics. Systems made of just two interacting atoms are 
especially relevant, being the building blocks of many- 
body strongly correlated states [3] and allowing for com- 
parison with fully understood theoretical models [H |5] . 

In a fascinating experiment [6 , the Heidelberg group 
has prepared in a one-dimensional trap (Fig. [l]) two ^Li 
atoms behaving as fermions of spin one-half and interact- 
ing through a tunable contact potential, g 5{xi—X2)' The 
Itt) (lii)) state has a fermionic-like orbital wave func- 
tion, being odd under particle exchange, '^^^(xi,X2) = 
— ?/;|l(x2, xi), whereas the W\) wave function is bosonic- 
like, ip^^{xi^X2) = '^^i{x2^xi). The time r spent by one 
atom to tunnel out of the trap is measured as a func- 
tion of the interaction strength g. As shown in Fig. |2j 



at ^ = oo the decay time r is the same for both W\) 
and Itt) (Itt) 7 being non-interacting, is independent from 
g). The coincidence of r is attributed to the Fermi-Bose 
duality (FBD, aka 'fermionization') [7]-[9], i.e., the en- 
ergy and wave function modulus of |tt) are the same as 

Itt), |V^n(^i'^2)|^ = |V^tt(^i'^2)|^. This duality was 
predicted long ago for the Tonks-Girardeau (TG) gas [7] 
and it has been demonstrated for repulsively interacting 
bosons [1^. At finite g both TG {g > 0) and super-TG 
{g < 0) states [3 [11] of |tt) are accessed, the former be- 
ing the ground state and the latter the first excited state, 
here stabilized by the absence of three-body collisions [6j . 

The experiment raises two issues: (i) A rigorous ar- 
gument for the coincidence of r was not given. In fact, 
the FBD applies only to those observables that involve 
the wave function square modulus [7|, whereas r origi- 
nates from the convolution of orbital states inside and 
outside the trap [12]. Indeed, in two dimensions r may 
depend on the wave function phase [131 El- (ii) So far 
r has been computed using a single-particle theory [6], 
including the effect of interactions through energy renor- 
malization [black dashed lines in Fig.|2|a)]. As the agree- 
ment with the experimental data (points) is qualitative, 
question arises whether the shape of the interacting wave 
function may affect r. 

In this Letter a theory of tunneling for two interacting 
atoms is developed, based on the calculation of the quasi- 
particle wave function (QPWF) [15l[16]. The QPWF is 
the interacting counterpart of the single-particle orbital 
resonating at the energy e of the atom tunneling out of 
the trap, dressed by the interaction with the companion 
atom. At ^ = cx), the QPWFs of both |tt) and |tt) 
are found to be identical in the tunnel barrier [Fig. [sj^f)], 
hence the decay times r coincide. At finite ^, r is eval- 
uated numerically from the exact wave function for har- 
monic confinement [4 plus perturbation theory (PT) to 
account for the anharmonic terms of the potential trap. 
The two-body correction significantly improves the agree- 
ment with experimental data [solid lines in Fig. [2|a)], 
illustrating a novel few-body effect, similar to orthogo- 
nality catastrophe [17]. 



2 



1000 



o 

CD 
C/D 

lr» 

CD 

E 

>> 
cc 
o 

CD 

Q 



3 



WKB 
QPWF 
WKB + PT 
QPWF + PT 



100 




1000 



100 



10 



1.9 



D) 

c 1 .5 

CD 
CD 

§ 1-1 

3 





1 1 1 1 


< 








7 (b) , ,1 


1 1 1 1 



80 

40. 



O 



-40 C 



-80 



720 750 780 810 840 870 900 
Magnetic field B (Gauss) 

FIG. 2. (Color online) (a) Decay time r vs magnetic field 
B. The points with error bars are the experimental data 
[6], the dashed and solid lines are respectively the WKB (ro) 
and QPWF predictions. The green [light gray] lines include 
the PT correction to the tunneling energy s. (b) Interaction 
strength g (red [gray] curve) and e (black and green [light 
gray] curve) vs B. g is taken from Ref. |6]and e is computed 
after Ref . |4j (black curve). The green (light gray) line includes 
the PT correction. The dashed lines are guides to the eye. 



The starting point of the theory is the Hamiltonian 



H = 



2m 



E 



dx 



2 + ^(^i) 



^g5{xi-X2), (1) 



with m being the mass and V{x) an effective one- 
dimensional potential. The plot of V{x) in Fig. [l] repro- 
duces the setup of Ref. |6| Atom-atom interaction is ne- 
glected outside the trap, that is approximately parabolic 
at low energy. Here the two atoms are either spin- 1/2 
fermions or spinless bosons. The treatment for bosons is 
the same as for hence only fermions are considered. 

Following Bardeen [12] , the space is divided into three 
regions: the trap (I), the barrier (II), and the vacuum 
(III) (Fig. [l]). The boundaries are the classical turnine; 
points Xa and X5, hence region II extends from Xa to X5. 
There are two classes of single-particle states: 
(i) Vacuum stationary waves Xei^) of (continuous) en- 
ergy £. These are plane waves in region III reflected 
at the barrier and dropping exponentially with distance 
into region II. For example, in the WKB approximation 



one has XeC^) = Ck{x) ^^'^ cos{j^k{x')dx' — ii / A) for 
X > X}) and X£(^) = Ck{x)~^^'^2~^ ex.'p {— j^\{x')dx') 
for Xa < X < X}), where C is a normalization constant 
and k{x) = [{2m/h^)\e-V{x)\]^/^. For x < Xa Xs{x) 
drops smoothly to zero, instead of oscillating, so Xei^) is 
not a good solution for x < Xa- 

(ii) Trap eigenstates 4)n{x) (n = 0, 1, 2, . . .), which have a 
vanishing tail in region II and smoothly go to zero in re- 
gion III. In region I (j)n{x) is approximately an eigenstate 
of the harmonic oscillator (HO) with energy ficJo(^+l/2). 
To summarize, trap states are solutions of the single- 
particle Schrodinger equation in regions I and II but not 
in III, whereas stationary waves are solutions in regions 
II and HI but not in I. Therefore, none of the above 
states are eigenstates in the whole space, but they are 
approximately orthogonal. 

Similarly, ^0(^17^2) and ^m,£(^i7^2) are states of the 
entire system with two fermions that differ in the transfer 
of an atom from region I to region HI. ^0(^17^2) has 
two atoms in the trap, either distinguishable (^0 = V^ti) 
or indistinguishable (^0 = V^tt or '^u)- The latter are 
non interacting as ?/^^^(x,x) = for symmetry. ^rn,e 
is the non-interacting state with one atom left in the 
trap orbital (j)m{x) and the other one transferred to the 
vacuum state Xei.^) (Fig. [l] depicts ^ce)- Thus ^0 is 
a solution of the interacting Schrodinger equation with 
energy Wq in region I and II but not in HI, and ^rn,£ is 
a solution with energy Wm.e in region II and HI but not 
in I. Both ^0 and ^rn,£ are good solutions in region II. 

Bardeen has shown [12 that the matrix element M^^e 
for the transition from ^0 to ^rn,e^ as long as the en- 
ergy is conserved {Wq ^ Wm,^), is Mm,e -h^ Jm.e/'^rn, 
with £ being proportional to the matrix element of 
the probability current density operator. 



Jm,s = / / dX2 S{Xi- 

•_i J — 00 J —00 



(2)' 

and Xbar being any point in the barrier region II. It is im- 
plicit in Eq. ([2| that ^0 and ^m,£ have the same spinorial 
component as tunneling does not affect spin. The decay 
rate 1/r may be estimated from Fermi golden rule. 



1 27r 



(3) 



The non-interacting case is easy to work out. The 
trap ground state ^0 is the Slater determinant of the 
(n = 0)(n = 1) configuration, 

V^n(^i'^2) = [^0(^1)01(^2) - 00(^2)^1(3:1)] . 

Among possible final states, the most relevant one is the 
in = {)){£ = 3hw/2) configuration ^o,^, 

^0,^(^1,^2) = [(l)o{xi)Xs{x2) - (l)o{x2)Xs{xi)] , 
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' (c) 

— ^Qp(^) 




FIG. 3. (Color online) Quasiparticle wave function (pQp{x) 
(black curves) and relative- mot ion wave function ^r(a:) (red 
[gray] curves) of |ti) at various g. (a) ^ = 0+. (b) 
g/(2^/^huJoiHo) = 3.14. (c) ^ = oo. (d) g/ = 
—3.77. (e) g = 0—. (f) The first excited state of the harmonic 
oscillator, 01 (x) (blue [dark gray] curve), being the QPWF of 
Itt)? is compared with (pQp{x) of at ^ = cxo. 



providing the transfer matrix element Jo,3hoj/2^ 



dx 



Xs{x) 



dx 



=3huj/2 



(4) 



Equation Q is the obvious non-interacting result for 
resonant tunneling between the highest occupied orbital 
of the trap, (/>i(x), and the vacuum state, X£=3huj/2{x)^ 
with the other spectator atom frozen in the orbital 0o(^) 
of the trap. In principle, there are final states other 
than ^o,3;iu;/2 allowed by energy conservation, like the 
(n = l){e = fux)/2) configuration. However, the corre- 
sponding matrix elements are negligible since wave func- 
tion tails drop exponentially with energy in the barrier. 

The case of two distinguishable fermions is not trivial, 
as ?/^^|(xi,X2) may not be written as a single Slater de- 
terminant. Nevertheless, the exact analytical solution is 
available in the relative- motion frame |Sl |5] . Explicitly, 

V^n(^l' ^2) = ^Cm([2:i + X2]/2)^r{xi - X2), (5) 

with ^CM and being functions of the center-of-mass 
and relative-motion coordinates, respectively. Since the 
singlet spin part is an odd function under particle ex- 
change, must be even, as shown in Fig. [3] (red [gray] 
curves), where ^r(^) is plotted for various interaction 
strengths g. Note that ^r(^) has two nodes for ^ < 
(super-TG regime). The wave function ^cm is just the 
lowest HO state in the center-of-mass frame. 

The final state ^o,£ is again a spin singlet, with an 
atom left in the trap lowest orbital (p^ix) and the other 
one transferred into the vacuum state Xel^)? 



^0,£(^l,^2) 



V2 



[(t)Q{xi)xE{x2) + (I)q{x2)Xs{xi)] . (6) 



Here the resonance energy e is fixed by energy conserva- 
tion, e is plotted in Fig. [2|b) (black curve) as a func- 
tion of the magnetic field B that controls the interaction 
strength g (red [gray] curve) in the experiment. 

By inserting expressions ([5| and (|6| into Q one ob- 
tains the important result 



(7) 



from 



with ipQ^{x) being the QPWF of the atom escaping 
the trap, defined as 

/oo 
dx'<t>o{x'Wcy,{[x' + x]/2)M/;(x' - x). 
-00 

(8) 

From the comparison of Eq. ^ and Q it is clear that, 
even in the presence of interaction, the tunneling is due 
to the resonance between the vacuum state Xei^) ^^^d 
the quasiparticle state ipQ-p^x). With respect to the non- 
interacting case of indistinguishable fermions, ipQ-pix) 
generalizes the role of the highest occupied orbital (j)i{x). 
Conversely, ^i(x) is the QPWF of the trivial case. 

At infinite interaction, ^ = 00, ^r(^) is 'fermionized' 
[red (gray) curve in Fig. |3jc)]. This means that ^r(^) 
overlaps with the modulus of the wave function of the 
non-interacting state |tt)- The latter is just |(/>i(a:)| in 
the relative- mot ion frame. The QPWF of W\) becomes 



2I/2 



^ ^HO 



V2^?io 



dx'e ^'^/^Ho(^' _ ^) + / dx'e ^'^/^ho (^ _ ^') 



,(9) 



with -^HO = {h/muJo)^^'^ being the HO length in the lab- 
oratory frame. From Fig. [sj^f) it is patent that (pQp{x) 
(black curve) is unrelated to (/)i(x) in the laboratory 
frame (blue [gray] curve), i.e., the QPWF of |tt). In- 
deed, (pQp{x) is an even function, with opposite peaks 
separated by a valley with significant weight in the cen- 
ter, whereas ^i(x) has a node. 

However, if x is in the potential barrier the first inte- 
gral between square brackets in Eq. (|9| may be neglected 
and the second one integrated over the whole space. This 
immediately provides (pQp(xbar) — ^i(^bar)- As it may 
be checked visually in Fig. |3jf), this approximate iden- 
tity is very well satisfied in region II. Since the energies 
e = 3J1W0/2 of both QPWFs are the same, it follows that 
the transfer matrix elements Q and ^ are identical and 
hence the decay times r, in agreement with the experi- 
mental data (circles with error bars) shown in Fig. [2|a). 

The inspection of Figs.[3]^a-e) over the whole range of g 
reveals that (pQp{x) is an amplitude distibution concep- 
tually different from ^^(x). This is evident from the very 
fact that the two functions have separate natural frames, 
respectively the laboratory and the relative-motion frame 
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(note the distinct lateral sizes in Fig. [3|. The form of 
^Qp(^) in the non interacting limits ^ = 0+ and ^ = 0— 
may be derived explicitly, providing respectively V2(j)o{x) 
[Fig. ga)] and M^)/V2 [Fig. (sfe)]. Overah, (^qp(x) 
smoothly interpolates between these two limits as 1/g 
goes from +oo to — oo. 

The magnitude of r depends on the tail of ipQp{x) in 
the potential barrier [Eq. Q]. Since there the QPWF 
profile is similar to those of non-interacting HO orbitals, 
it may seem that interaction affects r only by renormal- 
izing the tunneling energy £. This is incorrect, as the 
norm Aqp of the QPWF, Aqp = f dx\(pQp{x)f ^ varies 
strongly with the interaction strength. In fact, Aqp de- 
creases monotonously from the value 2 at 1/g = +oo 
down to 1/2 at 1/g = — oo, being one at 1/g = 0. 

Therefore, the decay rate 1/r of Eq. ([3| may be ap- 
proximately written as 



r To 

with the following provisos: (i) only one final state is 
considered, resonating at the energy e fixed by the inter- 
action strength g; (ii) tq is the decay time obtained in 
the absence of interaction for elastic tunneling at energy 
e; (iii) Aqp is the norm of the QPWF whose energy is s. 

Figure [2{a) shows the single-particle decay time tq 
(dashed lines) as a function of B (interaction strength) 
evaluated through the WKB formula used in Ref. [61 

In this approximation tq depends only on e and the shape 
of V{x). For Itt)^ ^0 is independent from B which was 
also observed in the experiment [red (gray) circles with 
error bars in Fig. [2|a)], being the state non-interacting. 
For Iti) (black dashed lines and circles) the agreement 
is only qualitative, with tq spanning three decades in 
the 1-1000 msec range and decreasing with increasing 
energy. The departure of tq from the measured time is 
systematic, with a blueshift on the repulsive side with 
the largest error of 45% at 5 = 725 G, and a redshift on 
the attractive side, with an error of 65% at 900 G. 

The QPWF correction based on Eq. ([To| improves 
significantly the agreement between the predicted [solid 
black lines in Fig. [2|a)] and measured (black circles) 
values of r on the positive side [JSI 1^ of the interac- 
tion strength ^, as the matching is almost within the 
error bars. On the negative side the estimated value 
of r is strongly blueshifted with respect to tq, with a 
residual error of 40% ai B = 900 G. Such discrep- 
ancy originates from the difference between the actual 
potential V{x) shown in Fig. [l] and the harmonic trap 
^Ho(^) = mujQ{x — xo)^/2 used to compute £ [4 , the 
higher the energy the stronger the anharmonicity. 

By employing PT to first order in V{x) — Vb_o{x)^ one 
may evaluate the correction to 6, as shown in Fig.[2|b) 
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FIG. 4. (Color online) Expansion of two-fermion states at 
various interaction strengths g on the basis of the Slater de- 
terminants obtained by filling with the atoms the orbitals of 
the harmonic oscillator (depicted as a ladder) and the sta- 
tionary waves in the vacuum. The numbers are the pertinent 
expansion coefficients. No label means a unity coefficient. 
The dashed red [gray] arrows point to the allowed tunneling 
transitions between initial and final configurations. 

(green [light gray] curve). Here the free parameter xq = 
0.659^HO has been fixed by minimizing As at ^ = oo 
and the resulting small value As = 0.02fia;o has been 
rigidly offset at any ^ < 0. Note that As <C hw^ in 
the considered range of 5, justifying the usage of PT a 
posteriori. The value of £ corrected by PT is then used to 
compute To through Eq. ( pTj ) [dashed green (light gray) 
line in Fig. [2|a)] and r through ([lO| [solid green (light 
gray) line]. Now the predicted value of r satisfactorily 
reproduces the experimental data, uncovering the two- 
body effect inherent in the loss of weight of the QPWF. 

The norm Aqp of the QPWF measures the matching 
of the initial two-atom state with that obtained by 
letting the current operator act on the final state ^o,£- 
To understand the evolution of Aqp with g it is useful 
to write in the laboratory frame. The simplest case 
of Itt) is schematized in the center column of Fig. [4] 
The initial state is a Slater determinant with one atom 
in 4)q{x) and the other one in ^i(x); the final state has 
one atom in (j)Q{x) and one in Xei^)- The coefficients of 
both determinants are one. The two states overlap if one 
transfers one atom from (j)i{x) into Xei^)^ indicated 
by the dashed arrow. This single-particle result may be 
taken as a reference, providing a unity contribution to 
the tunneling matrix element. This discussion holds also 
for W\) at ^ = oo. 

The continuous evolution of Aqp from ^ = — oo to 
g = {)— may be understood by considering the limit of 
^ = 0— , depicted in the right column of Fig. [I] (not acces- 
sible in the experiment). Three different configurations 
span ^0 in the lab frame: the two atoms either doubly 
occupy (/)i(x) or singly occupy (j)Q{x) and (j)2{x) with al- 
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ternate spins. However, only the latter configurations 
are effective in transferring the atom from ^2(^) into 
Xeix)^ providing two non-zero spin-conserving contribu- 
tions, each of weight 1/2 x 1/V2. The total contribution 
to the squared matrix element is 1/2. Such decrease of 
the QPWF weight, arising from the multiple configura- 
tions that generically span a many-body state, is familiar 
from the field of strongly correlated electrons interacting 
through Coulomb repulsive forces [15 . 

The positive branch of g shows a more exotic behav- 
ior. In the limit g = 0+ (left column of Fig. [4]) is 
a non-interacting configuration with 0o(^) being doubly 
occupied. This single configuration contributes twice to 
the final state, as both atoms may tunnel with opposite 
spms (2 X 1/^2). Therefore, the square matrix element 
is doubled. This coefficient may be regarded as a bosonic- 
like enhancement factor. Indeed, 7/^||(xi,X2) is the wave 
function of two spinless non-interacting bosons, both oc- 
cupying the level 0o(^)- 

In conclusion, the quasiparticle theory of tunneling 
here developed explains the available experimental data 
for two ^Li atoms. The comparison with measured decay 
times discloses a genuine two-body effect, related to the 
actual form of the interacting wave function. Therefore, 
such measurements may be regarded as a novel type of 
spectroscopy for cold atoms, paralleling successful tech- 
niques of condensed matter, such as photoemission and 
single-electron tunneling. 
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their experimental data prior to publication. I gratefully 
acknowledge support by Fondazione Cassa di Risparmio 
di Modena through the project COLDandFEW and 
from CINECA through CINECA-ISCRA project no. 
HP10BIFGH8. 
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